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ABSTRACT 

In 2009, Cygnus X-3 (Cyg X-3) became the first microquasar to be detected 
in the GeV 7 -ray regime, via the satellites Fermi and AGILE. The addition of 
this new band to the observational toolbox holds promise for building a more 
detailed understanding of the relativistic jets of this and other systems. We 
present a rich dataset of radio, hard and soft X-ray, and 7 -ray observations of 
Cyg X-3 made during a flaring episode in 2010 May. We detect a ~3-d softening 
and recovery of the X-ray emission, followed almost immediately by a ~1-Jy 
radio flare at 15 GHz, followed by a 4.3<r 7 -ray flare (E > 100 MeV) ~ 1.5 d later. 

The radio sampling is sparse, but we use archival data to argue that it is unlikely 
the 7 -ray flare was followed by any significant unobserved radio flares. In this 
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case, the sequencing of the observed events is difficult to explain in a model in 
which the 7-ray emission is due to inverse Compton scattering of the companion 
star’s radiation field. Our observations suggest that other mechanisms may also 
be responsible for 7-ray emission from Cyg X-3. 

Subject headings: black hole physics — X-rays: binaries — X-rays: individual 
(Cygnus X-3) 


1. Introduction 

The X-ray binary (XRB) system Cygnus X-3 (hereafter Cyg X-3), discovered by Giae- 
coni et al. (1967), is notable for its Wolf-Rayet (WR) companion (van Kerkwijk et al. 1996), 
short (4.8 h) orbital period (Bonnet-Bidaud & Chardin 1988), and dramatic radio variability 
(Waltman et al. 1994). The relativistic jets that it produces (Geldzahler et al. 1983) classify 
the source as a “microquasar” system. The distance to Cyg X-3 is ~7 or ~9 kpc, depending 
on the method used (Predehl et al. 2000; Ling et al. 2009). It is also uncertain whether the 
compact object is a black hole ( e.g . Cherepashchuk & Moffat 1994) or neutron star ( e.g . 
Stark k, Saia 2003). 

Cyg X-3 became more notable in 2009 when it became the first microquasar to be 
detected in the GeV 7-ray regime, via the satellites Fermi (Abdo et al. 2009) and AGILE 
(Tavani et al. 2009). A 4.8-h modulation in the Fermi data makes the detection definitive. 
Emission of 7-rays has been detected by both observatories during every “high/soft” X-ray 
state to occur since they began science operations (Abdo et al. 2009; Tavani et al. 2009; 
Corbel & Hays 2010; §3). The high/soft state is also associated with radio flaring and is 
closely tied to the launching of relativistic jets {e.g. Fender et al. 2004; Koljonen et al. 2010). 

Various physical processes are expected to be involved in the emission and propagation 
of high-energy photons in Cyg X-3. Inverse Compton (IC) upscattering of the intense (~10° 
erg cm -3 ) radiation field of the WR companion by high-energy (7 > 10 3 ) electrons in the 
relativistic jet is the most efficient source of GeV emission (Abdo et al. 2009; Dubus et al. 
2010). If the density of relativistic protons in the jet is sufficiently high, inelastic collisions 
between these and protons in the stellar wind will produce neutral pions that decay into 
7-rays (among other products) as well (Romero et al. 2003). In this case, there will be 
further 7-ray emission due to secondary leptons created during the hadronic interactions 
(Orellana et al. 2007). The typical hadronic content of microquasar jets is unknown: large 
numbers of hadrons impose significant energetic constraints, but there is direct evidence of 
their presence in at least one system, SS 433 (Kotani et al. 1994, 1996; Migliari et al. 2002). 
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A 7 -ray flare due to IC upscattering by jet electrons will be followed by a radio flare a few 
days later as those electrons cool, rarify, and become transparent to their radio synchrotron 
emission. This phenomenon has been observed repeatedly but not invariably (Abdo et al. 
2009; Tavani et al. 2009). If the 7 -ray emission is primarily hadronic, the result will be 
substantially the same, as secondary processes will create a similar population of energetic 
leptons. Abdo et al. (2009) determine an approximate 5 ± 7 d lag between radio and 7 -ray 
emission. Inspection of their Fig. 2 and the data presented by Tavani et al. (2009) suggests 
that this uncertainty is primarily systematic rather than statistical; i. e . , that the radio / 7 - ray 
lag varies significantly from event to event. Blazar observations often show lag variability, 
both between sources and in the same source over time, with theoretical support for sign 
changes in the lag between 7 -ray and other bands (Aller et al. 2010). 

In 2010 May, multiple monitoring programs detected that Cyg X-3 was becoming more 
active in 7 -rays and was entering the high/soft state (Bulgarelli et al. 2010; Corbel & Fermi 
Large Area Telescope Collaboration 2010; Kawai et al. 2010; Corbel & Hays 2010). We 
describe multiband observations of Cyg X-3 made during this episode (§2) and present our 
results (§3). Our data show a 7 -ray flare apparently lagging a radio flare, a sequencing 
opposite of that expected in the typical interpretation. We discuss and interpret this result 
(§4) and, finally, present our conclusions (§5). 


2. Observations & Data Analysis 

We observed Cyg X-3 during its 2010 May flaring episode with a variety of instruments. 
Below, we describe the observations and their analysis. Many dates and times we mention 
are near MJD 55,340. For brevity, we express these in an offset MJD, defining OJD = 
MJD — 55, 300. Thus OJD 0 is JD 2,455,300.5 or 2010 Apr 14.0 UT, and the 7 -ray flare was 
detected on OJD 43 = 2010 May 27 UT. 


2.1. Radio 

The Arcminute Microkelvin Imager (AMI) arrays (Cambridge, UK) are two aperture 
synthesis telescopes mainly used to study the cosmic microwave background (Zwart et al. 
2008). The observations described herein were made with the Large Array (AMI-LA), the 
reconfigured and reequipped Ryle Telescope, consisting of eight 13-m antennas with a max- 
imum baseline of about 120 m. observing in the band 12-18 GHz. The angular resolution is 
typically 25". Monitoring of small-diameter sources is undertaken as described in Pooley & 



Fender (1997); observations are interleaved with those of a phase-reference calibrator, and 
after calibration the data for individual baselines are vector-averaged. The in-phase compo- 
nent then provides an unbiased estimate of the target source’s flux density. The amplitude 
scale was calibrated by (at least) daily observations of 3C 48 and 3C 286, both of which are 
believed to be very nearly constant on long timescales. 

Should there be emission on a scale resolved by some of the baselines, the flux density 
estimate would be incorrect. There is an extended region of low-brightness emission around 
Cyg X-3 (Sanchez-Sutil et al. 2008) which is detected on the shortest baselines at extreme 
hour angles with the AMI-LA, and is sufficiently bright to need separate treatment when 
the central source in Cyg X-3 is “quenched” (flux density < 10 mJy at 15 GHz). No such 
correction was required in the observations described here. 

Observations of Cyg X-3 were also made with the RATAN-600 telescope (Nizhnij Arkhyz, 
Russia) as part of an ongoing microquasar monitoring campaign (Trushkin 2000; Trushkin 
et al. 2006). The observations and data analysis were performed as described in Trushkin 
et al. (2006). In this work we show results at 11.2 GHz; near-simultaneous observations were 
made at 4.8 and 7.7 GHz as well. 

Cyg X-3 was observed by the Allen Telescope Array (ATA; Northern California, USA; 
Welch et al. 2009) six times in the period OJD 26-43 as part of a larger transient search. 
Continuum images at 3.09 GHz were made after subtraction of model of the static sky. 
Source fluxes were determined in ^10-minute segments by fitting point source models in the 
image domain. 


2.2. X-Ray 

Our study includes hard X-ray observations made with the International Gamma-Ray 
Astrophysics Laboratory ( INTEGRAL ) satellite (Winkler et al. 2003). We used the INTE- 
GRAL Soft Gamma-Ray Imager (ISGRI; Lebrun et al. 2003) instrument and included two 
observations that were part of a program to observe Cyg X-l (PI: Wilms). Although the 
primary target of the observations was Cyg X-l, the ISGRI FOV is large enough to in- 
clude the entire Cygnus region, and our target was in the FOV throughout the observations. 
The first observation occurred during revolution 929, starting at OJD 40.625 and ending at 
OJD 41.940. The second observation occurred during revolution 938, starting at OJD 67.469 
and ending at OJD 68.805. A preliminary report of the first observation was given in Tomsick 
et al. (2010). 

We reduced the ISGRI data using the Off-line Scientific Analysis (OSA v9.0) software 
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package. The program ibis_science_analysis is the primary tool for extracting the stan- 
dard data products. We produced and inspected the ISGRI image in the 20-40 keV band, 
and then made energy spectra and light curves for all of the bright sources in the field. After 
correcting for instrumental deadtime, the total exposure times for Cyg X-3 from the first 
and second observations are 62,420 s and 66,540 s, respectively. We produced 20-40 keV 
light curves with a time resolution of 1 ks. 

We also obtained one-day average quick-look measurements from the All-Sky Monitor on 
the Rossi X-Ray Timing Explorer (RXTE/ ASM; Levine et al. 1996) and transient monitor re- 
sults from the Burst Alert Telescope on the Swift j - ray burst mission (Swift/ BAT; Barthelmy 
et al. 2005). These results are provided publicly by the RXTE / ASM and Swift / BAT teams, 
respectively. 


2.3. y-Ray 

Fermi Science Tools 9.17 and HEASoft 6.9 were used to reduce and analyze all Fermi- 
LAT (Atwood et al. 2009) observations within 10° of Cyg X-3 that took place in the range 
OJD -59-150 (2010 Feb 14 - 2010 Sep 11 UT), providing a baseline of ~100 d before and 
after the flare of OJD 43-44. Cyg X-3 is within the LAT FOV for ~15 ks per day. Given 
this and the typical flux of Cyg X-3 in the LAT band, we were compelled to use 1-day time 
bins in the following analysis. 

Thirty arcminutes away from Cyg X-3 is a comparatively bright pulsar, PSR J2032+4127. 
Using gtephem and gtpphase, we added the phase from the most recent ephemeris available 
(Abdo et al. 2009) to the events file. Then, we extracted events corresponding to the off- 
pulse phases (0-0.12, 0.2-0. 6, and 0.72-1; Camilo et al. 2009). This removed 20% of the live 
time, so the exposure time of the off-pulse events file was corrected accordingly. 

With gtselect and gtmktime, good time intervals from the off-pulse events file were 
selected in the 0.1-10 GeV energy range. At high energies (>10 GeV), Cyg X-3 is expected 
to have negligible emission. The events class was set to !i 3” which selects only high-quality 
diffuse-background photons. To minimize background albedo photons from the Earth’s limb, 
zenith and rocking angles were restricted to less than 105° and 52°, respectively. 

Exposure maps were generated by gtexpcube and gtexpmap while gtbin created photon 
counts maps in the region of interest. Emission from Cyg X-3 was not apparent on the counts 
maps, which is unsurprising given that Cyg X-3 is a faint source wdiose peak emission is at 
the low end of the LAT band where the point spread function is 5-10°. Complicating matters 
is the high level of diffuse background emission. 
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Following the procedure described in Abdo et al. (2009), unbinned likelihood analysis 
was performed for the daily bins with gtlike considering photons inside a 7° radius of 
Cyg X-3 from all 1FGL (Abdo et al. 2010) sources up to 5° away, and all bright (detection 
significance > la and flux [>100 MeV] > 5x 10 -8 ph cm -2 s -1 ) 1FGL sources up to 20° away. 
The instrument response function was “Pass 6 v3” (P6_V3_DIFFUSE), and the convergence 
relied on the NEWMINUIT method. Spectral models and input parameters from the 1FGL 
catalog were used except for pulsars 1FGL J2021. 0+3651 and 1FGL J2021. 5+4026 which 
had exponential cutoffs in addition to their power laws. Besides Cyg X-3, whose photon index 
and normalization were free, all sources had spectral parameters fixed to the 1FGL values. 
We included models for the galactic (gll_iem_v02) and extragalactic (isotropic_iem_v02) 
diffuse emission. The normalizations were left free to vary, as was the photon index of the 
galactic component. 


3. Results 

Long-term lightcurves of Cyg X-3 in soft and hard X-ray bands are presented in the 
top panel of Figure 1. The system was in a soft X-ray state during OJD ~20-55, with an 
episode of particular softness occurring during OJD 39-41. 

The Fermi source flux integrated over 0.1-10 GeV is presented in the lower panel of 
Figure 1. In bins where the test statistic (TS ~ —2 In L, where L is the ratio of the likelihood 
of models without and with the source) is less than one, we compute and plot a lcr upper 
limit. A peak is seen starting on OJD 43. The TS value in this bin is 18.4, which translates 
approximately to 4.3 a, and the flux is (4 + 1) x 10 -6 ph cm -2 s -1 . We are unable to achieve a 
more precise timing of the flare. Shifting the binning by half a day results in lower-significance 
and lower-flux detections in the range OJD 42.5-44.5. Analysis using higher-cadence binning 
does not yield significant detections of the peak. 

Radio measurements from just before the y-ray flare are presented in the top panel of 
Figure 2. A ~1 Jy radio flare occurred on OJD 42.25. There is evidence of an earlier flare on 
OJD 41.11, with a maximum observed flux density of 366 mJy, but only the rising portion 
of the flare is observed, so its properties are poorly-constrained. As is shown in Figure 3, 
the 15 GHz flux density of Cyg X-3 is depressed in the seven days following the y-ray flare, 
with a mean flux density of 48 mJy and an 80 th -percentile observed flux density of 66 mJy. 
(That is, 80% of the measurements during this period are < 66 mJy.) Starting on OJD 51, 
the typical radio flux density increases by a factor of ~2 to ~100 mJy. 


The bottom panel of Figure 2 shows high-cadence hard X-ray measurements from before 
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Fig. 1. — Long-term lightcurve of Cyg X-3 in X- and 7-rays. Shaded region: timespan 
depicted in Figure 2. Vertical black line: time of maximum observed 15 GHz flux density. 
On some days, our derived Fermi upper limits are smaller than the canonical background 
rate, 0.3 x 1CT 6 ph cm -2 s -1 ; on those days, we use the background rate as the upper limit. 

the 7-ray flare, made as Cyg X-3 exited the softest portion of its soft state. They show the 
recovery of the hard X-ray flux, subject to orbital modulation. The modulation makes it 
difficult to determine the significance of the drop in flux seen in the last three measurements; 
it may be evidence of a significant reduction in hard X-ray emission precursing the radio 
flare. Combining the Swift / BAT and INTEGRAL measurements, the hard X-ray flux has 
recovered by OJD ~42.0, about 5 h before the observed radio flare. The 7-ray flare then 
occurred 2. — 2 d later. 


4. Discussion 

If our observations are taken at face value — that is, one assumes no significant ra- 
dio activity during coverage gaps and that the flaring stems from one ejection event — the 
sequencing of the radio and 7-ray emission is inconsistent with the companion IC model 
discussed in §1. Emission in 7-rays following a radio flare could be interpreted as a reener- 





OJD = MJD - 55300 


Fig. 2. — High-cadence lightcurve of Cyg X-3 in radio and hard X-rays in the time immedi- 
ately before the 7-ray peak (— OJD 43-44). The X-rays are subject to orbital modulation 
with a 4.8-h periodicity. Vertical blue markers : times of X-ray minima according to parabolic 
ephemeris of Singh et al. (2002). Shaded region : partial timespan of Fermi peak. 

gization of the relativistic jet electrons by a shock (Dubus et al. 2010). With a typical jet 
speed of ~0.5 c (Tudose et al. 2007) and a delay of ~1 d, the reenergization would occur 
~100 AU from the system. This is much larger than the orbital separation but three orders 
of magnitude smaller than the distance at which the jet is expected to interact significantly 
with the ISM (~1 pc). Such reenergization could be due to a collision between the jet and 
a dense clump of the WR wind (Araudo et al. 2009). In this scenario, the absence of 7-ray 
emission at the time of ejection can be explained by absorption within the system (Cerutti 
et al. 2011), although this model must be reconciled with the observations of Tavani et al. 
(2009) and Abdo et al. (2009). The reenergization would lead to additional radio emis- 
sion as the 7-ray-emitting electrons cool, although the ambient magnetic field, and hence 
synchrotron luminosity, will be much weaker than that found close to the system. 

If the 7-ray flare is interpreted as the result of a discrete ejection event, the lack of 
a notable subsequent radio flare could be explained by the ejected material being largely 
hadronic. While hadronic 7-ray emission is not as efficient as IC upscattering, it does yield 
a higher ratio of 7-ray to radio luminosity (Romero & Vila 2008). The secondary leptons 




-9- 



OJD = MJD - 55300 


Fig. 3. — Radio lightcurve of Cyg X-3 during and after the 7-ray flare. The three panels 
proceed in chronological order from top to bottom. Red arrow, time of maximum observed 
15 GHz flux density (= 926 mJy). Shaded region: timespan of the Fermi peak (= (4 ± 1) x 
10~ 6 ph cm" 2 s _1 ). Dashed box: period of quiescent radio emission. 

due to hadronic interactions would, however, radiate, and detailed studies typically find that 
their bolometric luminosity is comparable to that of the primaries (e.g. Vila & Romero 2010). 

Should our data be taken at face value? Application of a simple synchrotron-cloud model 
(van der Laan 1966) supports the intuition that radio flares with sizes comparable to that 
of the largest observed could have occurred without detection, wflth the modeled lifetimes 
being ~0.2 d. While this model has had ambiguous success when applied to Cyg X-3 (Fender 
et al. 1997), its simplicity is advantageous for our sparse data, and it has been successfully 
applied to observations of other systems (e.g. Fender et al. 1997: Wilms et al. 2007; Prat et al. 
2010). To obtain a better understanding of the radio behavior of Cyg X-3 we compared the 
2010 May data to 15 years of archived AMI-LA / Ryle Telescope observations. The Jy 
radio flare is unusual: no flux densities > 400 mJy are detected in a 500 d span around the 
observations we describe. 
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In light of the predictions of the companion IC model, we consider particularly the 
radio observations after the observed 7-ray flare. The duty cycle of the observations is only 
6.8% in the seven-day post-flare quiescent interval. We identified 51 epochs similar to this 
in the archives, each lasting at least 7 days and having an 80 th -percentile observed flux less 
than 70 mJy (cf. §3). These archival measurements have a duty cycle comparable to that 
of the 2010 May data (4.9%), but a much larger time on-source (55.5 days). Because this 
time is large compared to the total duration of the 2010 May post-7-flare quiescence, a 
search for flares in the archival dataset can constrain the likelihood of there having been an 
unobserved flare in the 2010 May dataset, assuming no long-term evolution in the behavior 
of the quiescent state and stochastic flaring. 

The epochs were identified by exhaustively searching for seven-day segments meeting 
the aforementioned criteria, then lengthening these segments as much as possible without 
violating the 80 th -percentile constraint. Segments separated by less than three days were 
merged, in a few cases shortening the segment somewhat to preserve the statistical constraint. 
This method of construction does not bias against epochs containing rapid flares, no matter 
their size. Qualitatively, the lightcurve of the typical epoch starts high, drops to very low 
flux densities, and then becomes high again, possibly with rapid flares in the middle. 

The largest flux density seen in the selected archival data was 524 mJy, in the context 
of a single rapidly-evolving (risetime ~0.1 d) flare that would not have necessarily been 
detected in the 2010 May observing. (Here we exclude slow, large, epoch-terminating flares 
that would have been easily seen in the 2010 May data.) Approximating the rate of such 
flares as one per 55.5 days of observing, we derive an 11% chance that such a flare occurred in 
the 7 days after the gamma-ray peak without being observed. Ignoring the epoch beginnings 
and endings, which contain high flux values by construction, we find that the flux density of 
Cyg X-3 is larger than 250 mJy only 0.3% of the time during its quiescent state. 


5. Summary & Conclusions 

We have presented observations of a 2010 May Cyg X-3 flaring episode in the radio, 
soft and hard X-ray, and 7-ray bands. Our data show a sequence of three notable events: 
a particular softening and partial recovery of the X-ray emission, a rapid ~1 Jy radio flare, 
and a 7-ray flare. Interpretation of the data is made more challenging by the sparse sampling 
of the radio data and the low cadence of the high-energy observations as compared to the 
rapidity with which the radio emission can evolve. Nonetheless, the observations we do have 
and comparisons to historical data challenge the interpretation that the 7-ray emission is 
due to IC upscattering of the companion radiation field by high-energy electrons leaving the 
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system in a relativistic jet. because there is no evidence for the subsequent radio emission 
that one would expect to see from these electrons. 

While the companion IC model of 7 -ray emission from microquasars is clear and com- 
pelling, the lack of consistent radio/ 7 -ray timing lags call into question its completeness. 
There could be more than one mechanism responsible for the 7 -ray emission of Cyg X-3, 
becoming more or less relevant in different circumstances, or the reprocessing of the 7 -ray 
emission by effects such as pair production could be more significant than commonly as- 
sumed. Detections of or limits to very high-energy 7 -rays (>100 GeV; e.g. Aleksic et al. 
2010) or neutrinos from Cyg X-3 would aid in the understanding of the processes relevant 
to the emission in the Fermi band. 

The power of multiband monitoring of Cyg X-3 promises to increase significantly with 
the addition of 7 -ray data to the set of available observations. As Cyg X-3 inevitably pro- 
duces more flares, frequent observations at all wavelengths are important to build a detailed 
understanding of the launching and propagation of its relativistic jets. Radio observations 
of <2 h cadence with near-continuous coverage would be ideal for establishing an unambigu- 
ous relationship between radio and 7 -ray flaring. Based on the phenomenology we observe, 
intensive radio observations should be triggered at the ends of very soft X-ray states without 
waiting for the detection of 7 -ray flares. 
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